Load packages

suppressWarnings(suppressMessages({
library(Seurat)
library(magrittr)
library(imager)
library(EBImage)
library(STutility)
library(magrittr)
library(dplyr)
library(DT)
library(kableExtra)
}))

SVZ.Visium <- readRDS("SVZ.Visium/SVZ.Visium.rds")

Normalization

suppressWarnings(suppressMessages({SVZ.Visium <- SVZ.Visium %>% 
  SCTransform(verbose = F) %>%
  RunPCA(verbose = F) %>%
  RunUMAP(reduction = "pca", dims = 1:20,verbose=F)}))

Idents(SVZ.Visium)="labels"
SVZ.Visium <- RenameIdents(SVZ.Visium,"svze18"="E17.5","svzp2"="P2","svzp12"="P12","svzp22"="P22")
SVZ.Visium$labels <- Idents(SVZ.Visium)
SVZ.Visium@meta.data[["labels"]]=as.character(SVZ.Visium@meta.data[["labels"]])

Visualize only the spots

ST.FeaturePlot(SVZ.Visium, features = c("Sox2", "Pax6"), ncol = 2, grid.ncol = 2, pt.size = 4, show.sb = FALSE)

## Le chargement a nécessité le package : viridis

Visualize normalized gene expression in spots within sections

FeatureOverlay(SVZ.Visium, features = "Sox2",
sampleids = 1:4,
cols = c("light grey", "mistyrose", "darkred"),
pt.size = 3,
add.alpha = TRUE,
ncol = 4, show.sb = FALSE,
pt.alpha = 0.1,label.by = "labels")

Visualize normalized gene expression in spots with section

FeatureOverlay(SVZ.Visium, features = "Lrp2",
sampleids = 1:4,
cols = c("light grey", "mistyrose", "darkred"),
pt.size = 3,
add.alpha = TRUE,
ncol = 4, show.sb = FALSE,
pt.alpha = 0.1,label.by = "labels")

Other visualization parameters

FeatureOverlay(SVZ.Visium, features = "Sox2", 
              sampleids = 1:4,
              cols = c("lightgray", "mistyrose", "red", "darkred", "black"),
              pt.size = 1.5, 
              add.alpha = TRUE,
              ncol = 4,label.by = "labels")